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Recent lattice QCD calculations of the baryon spectrum are outlined. 

1 Introduction 

Quantum Chromodynamics (QCD) provides an excellent description of na- 
ture; however, the theory suffers from divergences that must be removed to 
render it finite. Lattice QCD provides an apriori non-perturbative regulariza- 
tion of QCD that makes it amenable to analytic and computational methods. 
No model assumptions other than QCD itself are needed to formulate the 
theory. This review surveys the rapidly evolving work in using Lattice QCD 
for calculations of baryon spectroscopy. Along the way, sources of systematic 
uncertainties in calculations are described and future directions are outlined. 

1.1 Regularization of QCD on a lattice 

As the starting point_for lattice QCD, the path integral formulation in Eu- 



continuum QCD are approximated with a discrete 4-dimcnsional lattice, with 
derivatives approximated by finite differences. Quarks are put on sites, glu- 
ons on links. Gluons are represented as 3 x 3 complex unitary matrices 
Un(x) = exp (igaA^x)) elements of the group SU(3) with vector potential 
A fi (x), coupling 5, and lattice spacing a. The vacuum expectation value of 
operators involves path integration over gauge and fermion fields 



The Gaussian integration over the anti-commuting fermion fields ip resulted 
in the det(M([/)) and M _1 (t7) factors with M(U) a lattice form of the Dirac 
operator. The gauge action Sq(U) approximates the Yang-Mills action of the 
continuum. The quenched approximation neglects the fermion determinant. 



clidean space 




The usual continuous space-time of 4-dimensional 
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The choice of working in Euclidean space resulted in no factors of i in 
the exponent multiplying the gauge and fermion actions. The path integral 
therefore resembles a 4-dimensional statistical mechanical model making it 
amenable to analytic as well as Monte Carlo methods for evaluation. 

Numerical predictions from lattice QCD are in principle exact (to some 
precision) after systematic errors are controlled. The statistical uncertainties 
go like 1/VN for N configurations of gauge fields in a Monte Carlo ensem- 
ble average. Systematic uncertainties include: (1) Finite volume - the lattice 
box must hold a hadron state, typically a lattice size of L ~ 2fm or more is 
needed. Several pion Compton wavelengths are needed m^L ~ 4. (2) Chiral 
extrapolations - calculations with small quark masses are expensive - extrap- 
olate observables to physical quark mass region (delicate!). (3) Discretization 
effects: inherent O(a) or 0(a 2 ) lattice uncertainty. One must extrapolate to 
continuum limit (a — > 0) to recover physical quantities. 



2 Confinement and Model Predictions - Static Quark 
Potentials 

A particularly useful role of lattice QCD is model testing. There is significant 
recent activity in the study of 3 quark potentials which provide phenomeno- 
logical insight into the forces inside a baryon. By gauge invariance, the quarks 
must be joined by 3 glue strings. These strings meet at a "gluon junction", 
which has been conjectured to be a non-perturbative excitation of the QCD 
vacuunJa. What is the area law behavior? One can test two ansatze. 

The F-ansatz predicts the potential grows linearly like V qqq cx o~ qq Ly 
where Ly is the minimal length of the 3 flux tubes necessary to join the .3 
quarks at the Steiner point. It is derived from strong coupling arguments^, 
and is consistent with the dual superconductivity confinement scenario. 

At large distances, the A-ansatz predicts instead that the potential grows 
linearly with the perimeter La of the quark triangle, e.g. V qqq oc a qq L^/2. It 
is derived_from a model of confinement by center vortices using a topological 
argument^. _ 

There is controversy as to which ansatztffl holds. Recent workQ claims 
that at short distances the potential approaches the A-ansatz but rises like 
the Y-ansatz at large distances. Departures from the A-ansatz appear above 
d qq ~ 0.7fm hence, the A model is more appropriate inside a hadron. However, 
recently Simonovu claims there is a field strength depletion near the F-junction 
which lowers the potential and could disguise the true behavior. Tests using 
adjoint sources could help reconcile the various claims. 



baryons '2002: submitted to World Scientific on February 1, 2008 2 



3 Hadron Spectrum 



3.1 Chiral Symmetry 

As mentioned before, for accurate lattice calculations systematic uncertainties 
need to be controlled. The discretization of the Dirac operator has been 
particularly troublesome since lattice QCD's inception and can significantly 
affect continuum and chiral extrapolations. The "doubling" problem is easily 
demonstrated by examining the lattice momentum representation of the free 
Dirac operator, namely X^7m^/j ~~ > » E M 7/' s ' n ( a ^)- The propagator has 
additional zeros at the momentum corners, e.g. ap M = 0, it so there are 16 
species of fermions in general. Originally, Wilson lifted the doublers by adding 
a Laplacian term that breaks chiral symmetry. In fact, the Nielson-Ninomia 
no-go theorems state one cannot avoid both doubling and chiral symmetry 
breaking with a local, hermitian action analytic in the gauge fields. This major 
theoretical problem has been solved with the recent advent of chiral fcrmion 
actions^ (e.g., Domain- Wall or Overlap fermions) and their use is crucial for 
matrix elements. How important is chiral symmetry for spectroscopy studies? 

Renormalization theory tells us that breaking a symmetry leads to in- 
duced quantum terms in an action. The Wilson fermion action has 0(a) 
scaling from the breaking of chiral symmetry. One can add a dimension 
5 operator (hyper- fine term) and rigorously improve scaling from O(a) to 
0{a 2 ). Scaling violations are dramatically reduced - mostly from improving 
chiral symmetry. Scaling violations are comparable with chiral fermion for- 
mulations. The conclusion is that chiral symmetry is important for accurate 
spectrum calculations^ at comparatively heavy quark masses. However, the 
benefits of chiral fermion actions with exact chiral symmetry are now being 
dramatically demonstrated as near physical quark masses are approached as 
will be shown below. 



3.2 Quenched Pathologies in Hadron Spectrum 

Clearly, precisely controlled lattice calculations come with the inclusion of 
the fermion determinant. However, because of their vastly reduced computa- 
tional cost quenched calculations are quite prevalent and one can gain impor- 
tant phenomenological insight into QCD, but the potentially large systematic 
errors induced in this approximation should be carefully ascertained. Sup- 
pressing the fermion determinant leads to well known pathologies as studied 
in chiral perturbation theoryld. There are missing vacuum contributions to 
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the disconnected piece of singlet correlators 

(ip(x)^ 5 ip(x)ip(y)^ 5 ip(y)) = (ir [75G (a;, y) 75 G (y, x)] c s ^ 

- N f (tt [75G (x, x) } C}S Tr [ 75 G (y, y)}^) . 

These effects are manifested in the 77' propagator missing vacuum contribu- 
tions with new double pole divergences arising of the form 

/ 7^4 eV ' X ( Tr 75G(x, x)Tr 75 G(0, 0)) = f P * 2 m 2 2 j 2 /„ + ... 
7 (27r) 4 + m£ p 771^. 

(1) 

How dramatic arc these quenched effects and to what extent do they affect 
the extraction of physical observables? One idea is to incorporate knowledge 
of quenched divergences in calculations and then attempt to extract useful 
information. 

3.3 Decay in the Quenched Approximation 

The m§ in Eq. ([!]) is the mass shift needed to recover the pseudoscalar sin- 
glet mass from the non-singlet pion. In LL3, the unique piece of the n' cor- 
relator - the hairpin - was computed directly. The parameter m§ was ex- 
tracted and one sees the lattice data is well described by the xPT prediction. 
With the shift, the rj mass (at non-zero lattice spacing) is determined to 
be 820(30)MeV with possibly large 0(a) scaling uncertaintities. A recent 
Domain Wall calculation^ gives 940(4)MeV. 

Further dramatic behavior is seen in the isotriplet scalar particle clq. 
There is an rj' — Tt intermediate state with missing contributions in the 
quenched approximation as shown in Fig. ([!]). In_£act, the oq correlator 
goes negative - a clear sign of violations of unitarity££l. From xPT, one can 
construct the aO correlator by including couplings between rj' — n states and 
rescattering states which can be resummed. The lightest aO correlator is fairly 
well described by a 1-loop resummed bubble term with n' mass insertion fixed, 
A mass m ao — 1.34(9)GeV was found. The new Domain Wall calculatioriia 
gives m ao — 1.04(7)GeV. The latter results does not exclude the possibility 
of a0(980) being a qq state. 

3.4 Quenched and Full QCD Hadron Spectrum 

The quenched low-lying hadron spectrum has been extensively studied by the 
CPPACS collaboration using Wilson fermionsEJ. Masses were computed at 
four lattice spacings and extrapolated to the continuum limit. Lattice sizes 
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Figure 1. Left: contributions to the ao propagator from an r\ — n intermediate state. Right: 
comparison of scalar ao propagator with the bubble sum formula fitted to the interval t=l-6 
(Ref£3). 

ranged up to 64 3 x 112 for a 3.2fm box. At each lattice spacing, the lightest 
pseudoscalar mass obtained was about 500 to 600MeV. Hadron masses were 
extrapolated in the quark mass via an ansatz motivated by the quenched X-PT 
prediction 

m 2 PS12 = A{m 1 + m 2 ){l - 5 [in (2Atoi/A^)] 

+ tn^/ — mi) In (mj/mi)} + B (mi + m2) 2 + O (to 3 ) (2) 
m H (mps) = mo + C 1/2 mp S + Cim 2 PS + C 3/2 m PS , C 1/2 oc 5 . 

The constants A, B, and Ci are fixed in the fits. Quenching effects were more 
clearly seen in the pseudoscalar channel. The calculation shows the basic 
hadron spectrum is well determined even in the quenched approximation to 
within 10% accuracy. The computational cost was roughly 50 Gigaflop-years. 

A subsequent two-flavor jdynamical calculation was made with four quark 
masses at 3 lattice spacingsti Box sizes range up to about 2.5fm. In the 
meson sector, the results are consistent with the original quenched calcula- 
tions and now agree to within 1% of experiment. Systematic deviations of the 
quenched calculation from experiment are seen demonstrating that sea quark 
effects arc important. The vector meson masses are increased after unquench- 
ing. This increased hyperfine splitting is consistent with the qualitative view 
that the spin-spin coupling in quenched QCD is suppressed compared to full 
QCD due to a faster running of the coupling constant. 

In the baryon sector, two-flavor dynamical sea quark effects are not as ap- 
parent. As seen in Fig. (^). The N and A masses are higher than experiment, 
but other masses are consistent. With only a 2.5fm box, finite- volume effects 
could well be large. Another concern is that the octet and decuplet chiral 
extrapolations have many parameters resulting in possibly underestimated er- 
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Figure 2. Baryon masses in two-flavor (filled symbols) mid quenched (open symbols) QCD. 
Both graphs have the same lattice spacing scale. (Ref.ta). 



rors and will be discussed more next. However, this is a significant calculation 
involving roughly one Teraflop-year of computations and well demonstrates 
the efficacy of lattice methods. 

3.5 Improved Chiral Extrapolations 

The Adelaide group-has been extensively studying higher order xPT effects on 
hadronic quantitiesli3. The basic upshot is that the naive chiral extrapolations 
in use are just too naive! In particular, they incorporate leading non-analytic 
behavior from heavy baryon xPT arising from B — » B'-k — > B intermediate 
states with B = N, A 

M B = a B + pBml + E s (m,, A) (3) 

where is a self-energy term. The coefficient of the m^. term is actually 
known analytically in contrast to Eq. (|2|). The basic argument is since xPT 
has a zero radius of convergence (or certainly not a well defined radius), a 
simple leading order approximation to Eq. ([}]) is quite a bad approximation 
at moderate m„. They use a simple regularization of the self-energy term. 

Fig. (^|) shows a comparison of a recent quenched and 2 + 1 dynamical cal- 
culation of the low lying hadron spectciim for a variant of staggered fcrmions 
at the same (physical) lattice spacing!!!! The various lines are from fits using 
Eq. (Q). One can see that the self energy term becomes significant for small 
lattice quark masses. The intriguing result is that the fit parameters cxb and 
Pb agree very well between the quenched and dynamical calculations. This 
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Figure 3. Left: contributions from various intermediate states to the quenched and un- 
quenched self-energy term E. Right: fit (open squares) to lattice data - quenched (open 
A) and dynamical (filled A) with adjusted self-energy expressions accounting for finite 
volume and lattice spacing artifacts. The continuum limit of quenched (dashed lines) rewnd 
dynamical (solid lines) are shown. The lower curves are for N and upper for A. (Ref. E3-) 



result can be used to justify the claim that the dominant effects of quenching 
is attributed to first order meson loop corrections. 

While quite intriguing it is fair to say there is some controversy over 
these results. At issue is the concern that once one uses any model to directly 
interpret lattice results, one has lost predictably. However, in defense once 
one used a chiral extrapolation at all one has chosen a model. Ultimately, 
the Adelaide's group work has shown that their is interesting structure in the 
"pion cloud" around a hadron and going to light quark masses is essential. 

3.6 Excited Baryons 

Understanding the N* spectrum gives vital clues about the dynamics of QCD 
and hadronic physics. Some open mysteries are what is the nature of the 
Roper resonance? Why is the ordering of the lowest-lying states - the positive 
and negative parity states - inverted between the N, A and A channels? 

The history of lattice studies of excited baryons is quite brief. Recently, 
new calculations are starting to appear using improved gauge and fermion 
actions. The nucleon channel is the most studied and work has focussed on 
two independent local interpolating fields 

Ni = Eijk (ufCj 5 dj) u k , N 2 = e ijk (ufCdj) j 5 u k . (4) 

Both interpolating fields couple to positive and negative parity states, so in 
practice parity projection techniques are used. Making the lattice anisotropic 
with finer discretization in time allows the behavior of the correlators to be 
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Figure 4. Left: effective masses for correlators with fields Eq. (^) corresponding to the 
nucleon (circles), its parity partner (diamonds) and tentatively the lowest positive parity 
excitation (cross) with an anisotropic clover action_( Rcf.Lj) • Right: masses in physical units 
obtained with the isotropic Overlap action (RefJiS). Solid symbols denote N(^ + ) states: 

ground (•) and lst-excited (*). Empty symbols denote JV(^ ) states: lowest (A) and 2nd 
lowest ([— ]). The experimental points (*) are taken from PDG. 



examined over many more time slices than on isotropic lattices. Additional 
tuning of the fermion action is needed to recover hypercubic symmetry. 

The left side of Fig. (|j) shows the effective mass for the nucleon and its 
parity partner on an anisotropic lattice using the clover actionLJ - an action 
improved to have reduced discretization uncertainties. Long plateaus (clean 
extraction of a mass) are seen demonstrating the efficacy of the method. For 
this quark mass, there appears to be the expected ordering of the states 
jy*i/2+ > jyi/2- > jyi/2+_ However, there is difficulty in approaching small 
quark masses, and the mass extracted with the N2 operator appears too large. 

A recent quenched calculatiorJIS using Overlap fermions (a chiral fermion 
action) going to much smaller quark masses reveals a dramatic decrease in 
the nucleon masses extracted only with the N± interpolating field as seen in 
the right side of Fig. (|4|). The authors claim the apparent crossing of the first 
excited N 1 / 2 and lowest TV 1 / 2- states is the demonstration of the physically 
correct ordering of states. However, it is also possible that at such light 
quark masses a decay threshold has been crossed and the mass observed is 
affected by missing dynamical effects via the mechanisms described in Sec. 3.3. 
If so, a finite volume check via shrinking the lattice box can reveal this. It 
is clear lattice calculations are really beginning to probe interesting excited 
state phenomena. With judicious use of finite volume techniques, physically 
relevant mass information can be extracted. 
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4 Conclusions 



First generation lattice calculations of excited baryon spectroscopy are ap- 
pearing. State of the art calculations require roughly 100 Gigaflop-year in 
quenched QCD and roughly 1 to 10 Teraflop-years in full QCD. The required 
resources are not available to the US lattice community The Dept. of En- 
ergy's SciDAC program is addressing this shortcoming and a large effort is 
ongoing in the U.S. to meet future computational needs. 
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